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We derive new uncertainty relation for one-dimensional mixed state given its purity and degree 
of non-Gaussianity. This relation extends the purity-bounded uncertainty relation for mixed states 
derived by V. V. Dodonov and V. I. Man'ko. For the special case of pure states it provides us with 
an extended version of the Robertson-Schrodinger uncertainty relation, saturated by a set of states 
which includes all the eigenstates of the quantum harmonic oscillator. We represent our results as a 
bound in a three-dimensional parametric space of mixed states and identify the regions of the bound 
realized by states with non-negative Wigner function. This takes us closer to a proper extension of 
Hudson's theorem to mixed states and permits us to visualize and compare the set of states with 
non-negative Wigner function and the set of states which minimize the newly derived uncertainty 
relation. 



I. INTRODUCTION 

The only pure states which saturate the Schrodinger- 
Robertson [l| uncertainty relation for the canonically 
conjugated coordinate and momenta, are the Gaussian 
states. The latter, due to Hudson's theorem 0, are also 
the only pure states with strictly positive Wigner func- 
tion Q. One may naturally wonder how these theorems 
can be extended to mixed states, and whether the "Gaus- 
sian" link between them persists. 

In the space of mixed states, different forms of un- 
certainty relation have been suggested varying with the 
choice of the quantities characterizing the mixed states 
i.e. purity, purities of higher order or various entropies 
(see [1] for a review). Among these proposals, a basic 
one, i.e. expressed in a closed form, is the purity-hounded 
uncertainty relation suggested by Dodonov and Man'ko 
. According to this relation the minimum uncertainty 
is a function of the purity characterizing the degree of the 
mixedeness of a state. In a more recent work Dodonov 
Q has proven that the mixed states which minimize this 
relation possess positive Wigner functions and that they 
are not extremely "far" from Gaussian states of the same 
purity. If instead of the purity, the von Neumann en- 
tropy is employed as a measure degree of mixedeness, 
the results are not identical. Bastiaans @ has derived an 
uncertainty relation based on the von Neumann entropy 
of mixed states, the entropy-hounded uncertainty relation 
and has proven that the mixed states of minimum uncer- 
tainty are the thermal (Gaussian) states. 

These results already give a preliminary answer to the 
questions we posed suggesting that, the link between 
the ability to minimize relevant uncertainty relations and 
positivity of Wigner functions seem to persist in the space 
of mixed states, though the Gaussianity is no longer a 
necessary requirement for extremality of the uncertainty. 
However, to be able to answer our questions in more com- 
plete way we need to have a closer view on the set of 
mixed states with positive Wigner functions. 

In a recent work [7| we have formulated an extension 
of Hudson's theorem, by identifying the maximum degree 



of deviation from a Gaussian state which can be reached 
by a mixed state with positive Wigner function, given its 
purity and its uncertainty. This problem does not have 
a simple solution and in we have only identified the 
bounds for any continuous classical distribution which 
may not necessarily correspond to a valid Wigner func- 
tion, and therefore the bounds are not tight for the set 
of quantum states. 

In view of this discussion, we extend in this work the 
purity-bounded uncertainty relation by adding one more 
parameter, the non-Gaussianity characterizing the dis- 
tance of the state from a Gaussian state with the same 
covariance matrix. In this way we are able to draw more 
complete conclusions on the overlap between states mini- 
mizing this uncertainty relation and the set of states with 
positive Wigner function. In addition, this extension 
permits us to show explicitly how the Gaussian states 
smoothly become extremal pure states in the context of 
uncertainty and positivity of Wigner's functions. 

In Secini we introduce the quantities which we choose 
to characterize one-dimensional mixed states and shortly 
review the already acquired results concerning uncer- 
tainty relations. In Sec. IIIII we derive new bound on the 
uncertainty of a mixed state characterized by its purity 
and non-Gaussianity. We name this relation the non- 
Gaussianity hounded uncertainty relation. In Sec. IIVI we 
identify numerically the parts of the bound realized by 
states with positive Wigner function, thus visualizing the 
overlap of the set of states with positive Wigner function 
and the states minimizing the non-Gaussianity bounded 
uncertainty relation. 

In a recent work [§j we have extended the entropy- 
bounded uncertainty relation using similar methods as 
we use here. However, that extension, as well as, the 
entropy-bounded uncertainty relation itself requires nu- 
merical methods to be expressed in a proper form. In 
Sec|Vl we discuss the main outcomes of this work and we 
compare them with those in [81] . 
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II. OUR PARAMETRIC SPACE 

For our purposes we need first to choose a convenient 
parametrization of the space of mixed states that wih 
permit us to formulate the extensions of the uncertainty 
relation and Hudson's theorem in a unified manner. 

Let us start with a general density matrix p deter- 
mining a quantum state. Apparently a measure of the 
degree of mixedeness of the state is an indispensable pa- 
rameter characterizing a mixed state. We choose as the 
most convenient measure for us, the purity ^ = Tr (p^). 

The second basic characteristic of the state p is its 
covariance matrix 7 defined as 

j,^^Tr{{{h~d,),{r,-d,)}p) (1) 

where f is the vector of the canonically conjugated ob- 
servables position and momentum (or quadrature opera- 
tors in quantum optics) f — {x,p)'^ , d =Tr{rp) is the dis- 
placement vector, and {•,•} is the anticommutator. We 
can put the displacement vector to zero with no loss of 
generality since the purity (as well as the quantities we 
will be interested in) is invariant with respect to d. 

The importance of the covariance matrix stems from 
the fact that it is measurable experimentally and that it 
is tightly connected, as we show below, to the uncertainty 
relation. Throughout the text we will consider quantum 
systems with only one pair of conjugate variables, like 
the system of one particle or one-mode state in quantum 
optics, therefore we consider 2x2 covariance matrices. 
The covariance matrix of p uniquely determines a Gaus- 
sian state pq which henceforth we will refer to as the 
reference Gaussian state for the state p and we choose as 
second parameter the purity of pa, 



fJ-G = i/y 711722 - I712I ■ (2) 

As a third parameter we introduce a measure of the 
non-Gaussianity, i.e., the distance of the state from its 
reference Gaussian state. At the first step we are go- 
ing to work with the trace overlap Tr {ppc) between the 
states p and pc ■ Although the trace overlap is not a mea- 
sure of distance, it can be employed together with pc: M 
for evaluating a correct measure of non-Gaussianity, the 
normalized Hilbert-Schmidt distance 01,0, 

c 1 rp . ^2^ AiG + M-2Tr(/5pG) 

^ = Tp^' ~Po))- ■ (3) 

Therefore, while a state p requires an infinite set of 
parameters to be fully described, it is enough for our 
purposes to work in the reduced parametric space {p, pa, 
S}. In [3] we have used the same parametric space. There 
we have formulated an extension of Hudson's theorem to 
mixed states by identifying an upper bound, though not 
tight, on the non-Gaussianity S for states with positive 
Wigner function. This bound is represented by a surface 
in the space {p, pa, 5}- 



Let us we note that the existing purity-bounded uncer- 
tainty relation Q can be represented as a line-bound in 
the plane {p,pg}- Before we justify this statement, we 
shortly review the historical evolution of the knowledge 
concerning uncertainty relations and for more details we 
address an interested reader to For pure states of 
covariance matrix 7, the first uncertainty relation 

711722 ~ 4 (4) 

(where h — 1) was proposed by Heisenberg [l^ and 
proven in its more strict form, 711722 ^ 1, by Kennard 
in 1927 . This uncertainty relation is saturated by 
coherent states and squeezed states whose major and mi- 
nor axis in the phase-space, are parallel to the x- and p- 
axis. Few years later, the Heisenberg uncertainty relation 
was generalized by Schrodinger and Robertson, indepen- 
dently, to a relation valid for any set of non-commuting 
observables and which, in the case of {x,p), takes the 
form 

\ll11l22 - \ll2\^ > 1- (5) 

The Scrodinger-Robertson uncertainty relation accounts 
for the possible correlations between the coordinates and 
therefore is saturated by all Gaussian pure states. Even 
though, the inequality Eq. ([5]) is valid also for mixed 
states it cannot be saturated by any mixed state. There- 
fore, for the sets of states with purity less than one, this 
uncertainty relation is not tight. A tight uncertainty re- 
lation for mixed states is the purity-bounded uncertainty 
relation derived by Dodonov and Man'ko 

\/ 711722 - |7i2|^ > "J" (a*) ■ (6) 

For mixed states the lower limit on the uncertainties de- 
pends on the purity of the state via a monotonic function 
$ satisfying 

$(1) = 1 and 1 < $(/i) < l/^x. (7) 

In view of Eq. ^ it is then evident that the purity- 
bounded uncertainty relation Eq. (jG)) is equivalent to a 
bound given by a line in the {p, pa} plane. Dodonov 
and Man'ko have proven Q that this bound is realized 
by states with phase-independent and positive Wigner 
functions and they have derived an exact but piece- 
wise expression for the function $ in Eq.®. In A we 
present an alternative, though parametric, expression for 
this line-bound whose construction permits for a gener- 
alization to the three dimensional parametric space. We 
also note here a quite unexpected fact visualized by the 
purity-bounded uncertainty relation. It is known that on 
one hand the pure states which saturate the Scrodinger- 
Robertson uncertainty relation for pure states are the 
Gaussian states. On the other hand the Gaussian mixed 
states (thermal states) maximize the von Neumann en- 
tropy 6] being, in a sense, maximally "disordered" (l2j 
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given Eq.(l2]). However, contrary to our intuition, 
the purity-bounded uncertainty relation dictates that the 
mixed Gaussian states are not the states which maximize 
the purity. 

In the next section, we extend the purity-bounded 
uncertainty relation by including into it the non- 
Gaussianity 5 of states 



711722 -I712I >F{^i,S). 



(8) 



For our derivation and in view of Eq. ^ it is convenient 
(as we explain in more detail in Sec lIIII) to conceive the 
inequality given by Eq.® as a bound in the introduced 
parametric space {fj,, iiqt S}. This bound valid for all 
mixed states is then compared, in Sec lIVi with the bound 
for states with positive Wigner functions. 



III. NON-GAUSSIANITY BOUNDED 
UNCERTAINTY RELATION 

In order to find function F in Eq.® it is enough to 
find the bound for the whole set of mixed quantum states 
in the parametric space {fi, fj,G,S}. It is obvious that 
such a bound will be a surface (since the space is three- 
dimensional) and if it is single- valued in the parameters /i 
and S it would represent the function F (/i, S). Our results 
show that the bound is in general single-valued, how- 
ever there is some particular domain of {//, d}. where the 
bound is double- valued. This problem does not emerge if 
one works with Tr (ppc) instead of S (see Eq. ([3])). That 
is why we derive our results employing the trace over- 
lap. With the help of Eq. ^ the results can be directly 
translated to the initial parametric space {/i, /xg,(5} and 
we do so, mainly for the visualization of the results. 

A possible way to proceed with the derivation of the 
bound is to search for the states which minimize or max- 
imize the purity fi if the values of the other two param- 
eters, HG and Tt{ppg), are kept constant. We conceive 
then the total bound consisting of two regions, namely, 
the region I where the purity p is minimized and the 
region II where the purity is maximized given pc and 
Tr (ppq). At the final step we recombine both regions and 
present the whole bound which we name non-Gaussianity 
bounded uncertainty relation as the abbreviation of the 
more exact title non-Gaussianity and purity bounded un- 
certainty relation. 



A. Region I: States of minimum purity 

Let us consider a general density matrix p character- 
ized by the purity pc of the reference Gaussian state pc 
and the trace overlap between p and pa, Tr (ppc)- By us- 
ing symplectic transformations we can set the covariance 
matrix 7 of p into a symmetric form where 711 = 722 and 
712 = 0. Symplectic transformations, such as squeez- 
ing and rotation, do not change the quantities of interest 



(see ref.Q) and therefore without loss of generality, we 
can assume that the reference Gaussian state for p, pc, 
is a thermal state. In the Wigner representation such a 

state, pg, is 



TT 



(9) 



where pQ is its purity defined in Eq. ([2]). 

The next step is to prove the Lemma stating that for 
a given state p with these characteristics, there is always 
a state of the same characteristics possessing a phase- 
independent Wigner function, that is of equal or lower 
purity than p. To prove the Lemma it is more convenient 
to work in the Wigner phase-space representation. In 
the general case, quantum state p possesses an angular- 
dependent Wigner function W (r, (p) and its purity can 
be expressed as 

P^2tt W{r,ipfrdrdip. (10) 



The trace overlap between p and pc is 

Tr (ppg) = 2tt W{r, ip) Wq (r) rdrdip. (11) 



Now let us construct a new state ps with phase-invariant 
Wigner function Ws (r) , by phase-averaging the Wigner 
function W (r, if), 



Ws (r) 



1 

2^ 



W (r, if) dip. 



(12) 



The reference Gaussian state for ps will be the same as 
for p, since phase-averaging cannot affect the angular- 
independent Wigner function Eq.®. The trace overlap 
among the symmetrized state p^ and pc will remain the 
same, as well. This is a straightforward result of substi- 
tution of the phase-independent Wigner function given 
by Eq. (IT2|) into Eq. (fTTj) . However, the purity ps of the 
symmetrized state ps is constrained to be smaller than, 
or equal to, that of p. Indeed, by applying the Cauchy- 
Schwarz inequality we have 

/is = 27r / / Ws (r)^ rdrdp 



W (r, p) W (r, <i>) rdrdpd^ 




W (r, (p) rdrdpd^ 



W (r, p) rdrdp = p. 




W (r, (f>) rdrdpd^ 



This concludes the proof of our Lemma. 

From this Lemma is straightforward to deduce that the 
bound corresponding to the states of minimum purity 
is achieved by states with angular-independent Wigner 
function and we proceed by applying the method of La- 
grange multipliers to identify these states. We note also 
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here that this Lemma is well in accordance with the re- 
sult of Dodonov and Man'ko (51 who have shown that for 
any given fic one can always choose a coordinate sys- 
tem where the states with angular-independent Wigner 
functions minimize the purity fi. 

Any state possessing an ang ular- indep endent Wigner 
function can be expressed jl3| as a convex combination 
of the number (Fock) states 



n=0 



n in 



(13) 



In order to identify the weights p„ of the states which 
minimize the purity 



71=0 



(14) 



we apply the method of Lagrange multipliers with the 
following constraints imposed on the solution, 



1. Normalization 



n=0 



(15) 



2. Fixed purity for the reference Gaussian state pc 
(an un-shifted thermal state as in Eq.® ) 



oo 

E 

n=0 



Pn (2n + 1) = l/^J.G■ 



3. Fixed overlap between p and pa 

2p,G (1 - Pg)" 



Tr {ppc 



n=0 



(1 + Mg) 



(16) 



(17) 



Applying the Lagrange multipliers method on Eg. p^ 
with the constrains EQS.([T5 |) - p7)) . we obtain the weights 
which realize the extremum solution 




and only if / (x) is convex and positive function in an 
interval [a;i,X2] of the positive semi-axis. When this con- 
dition is satisfied, the integer parts of the roots xi, X2 of 
the equation / (x) = define respectively the lower and 
upper limits of summation in Eq. (1131) . In the case that 
the equation has only one positive root we denote this 
root as X2 and put xi equal to zero. With this ansatz, 
the density operator that maximizes p in Eq. (jl4p under 
the constraints Eqs. ((T5|) - pT|) . can be re-expressed as 



P = 



(20) 

with rimin — rimax = [2^2] • (whcrc [x] gives the 

smallest integer that is greater or equal to x, and [x\ 
gives the integer part of x). The rank of the solution 
rank (p) — fimax — J^min .This ansatz can alternatively be 
interpreted as adding two unknowns, nmin and nmax, into 
the initial extremization problem and setting two more 
constrains, f{xi) = and f{x2) = 0, with f{x) given by 
Eq. ( [in]) , but for brevity reasons we do not repeat the 
procedure. 

Now, we are ready to determine the coefficients Ai. 
For every positive integer value of nmi„ (including zero) 
we express the unknown Lagrange multipliers Ai, A2 and 
A3, in terms of pc and X2. More precisely we employ the 
conditions Eqs. (fT5|) - (|T6| . and (IT9l) which look now as 



iAi+ A2n + A. 



PG (1 ~ P gT 

(I + A^g) 



n+1 



(21) 



^ iA,+A2n + As^ 



PG (1 - PgT\ ,r, , ^^ 1 / 

- (2n + 1) = 1/pG 



+ A^g)' 



Ai + A2X2 + A: 



PG (1 - paf 



(22) 



(23) 



p,^ = A,+A2n + A,^P^J^ (18) 
(1 + Mg) 

with the coefficients Ai determined by the conditions 1 — 
3. 

In order to restrict ourselves to positive density ma- 
trices we replace the index n in Eq. ([T8l) by a continuous 
variable x and introducing instead of p„ the function 

f{x)=A,+A2x + A,(^^^^^^. (19) 

The same ansatz is employed in Appendix A to derive 
an alternative expression for the purity-bounded uncer- 
tainty relation Eqs. (|5^ - (I55|) . One can see that the so- 
lution Eq. (jl8l) . corresponds to a valid density matrix if 



This system of equations in a more explicit form (after 
the summation over n) is presented in Appendix B. For 
each non-negative integer rimin one has to solve this lin- 
ear system of equations. After this, the last step is to 
substitute the derived expressions for Ai into Egs. p^ 
and (|17p and obtain parametric expressions for extremal 
/i^^ and Tr(pG'P)'^^ via pG and X2- In order to assure that 
JT-min is the same in all three equations, one has to require 
that the constrains / (rimin) > and / (rimin — 1) < are 
satisfied. 

The exact solution presented in the Appendix B is 
rather cumbersome. However, one can derive an approx- 
imate solution simply by setting [X2J = X2- It turns out 
that for riniin = one arrives to the following parametric 
expression 
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^ (8(20:2 
^ ((2: 



1) mg(mg + {'2x2^lG + A^G - 3) + y^"^ - {2x2 + 1) Vg^ + 4 (4a;i + Gxj - l) mg') 
^2 ((2a:2 + 2x2 ~ 1) Mg^ + (4a;2 + 2) /ic + 3) + (;Ug + 1) (2x2Mg + Mg - 3)) 2 
/iG (18 {2x1 + 2a;2 + l) Mg' + 12 (2x2 + 1) mg " 9 + (a^g + 1)' (2a;2MG + Mg - 3) 
^ (y-2((2a:2 + 2x2-l)AiG2 + (4x2+2)^G + 3) + (AtG + l)(22;2MG+AiG-3))2 ^ ^ 

. .ex ^ A^G (-4 (2x2 + 1) MG?/"^ + (a'g - 1)2/'"^ (2a:2MG + Mg + 3) - (mg+1) (2x2/^g + Mg - 3)) 
''^^^''^ y-^((2x2 + 2x2-l)AiG2 + (4.X2 + 2)MG + 3) + (/iG + l)(2x2AiG + AiG-3) ^ ' 



where y = (1 + AtG)/(l — Mg), 2:2 > 2. These equations 
represent the global solution in very good approximation. 

In Fig. [T] we present the derived solution which we 
caU region I of the bound, projected on the planes, (a) 
{AiG,Tr (ppg)}, and (h){nG-,^]- In Fig. [DJa), the whole 
bound is confined between the dotted straight line which 
corresponds to the Gaussian states and the solid line 
which represents the solution of rank 2 realized by states 
that are mixtures of any two successive number states 

Pr2 = a \n) {n\ + (1 - a) |n + 1) {n + l\ . (26) 

In Figlljb), the lowest dashed line represents the purity- 
bounded [5| uncertainty relation. One can see it is the 
lower limit of the more general bound derived here. The 
upper limit of the so far derived bound (the upper solid 
line) is the rank 2 solution Eg. ((26)) . This curve reaches 
the line ^ = 1 oi pure states (not shown in the graph) 
only at the points realized by the number states (indi- 
cated on the graph as n = 1,2,...). The area between 
this line and the /i = 1 line remains not covered by the 
derived bound. 

In the same figure (Fig[ljb)) one can notice that solid 
lines form loops which are arranged in horizontal rows. 
These lines define the borders between different ranks 
rank {p) = rimax — 'T-min of the solution in Eg. ([20]) . Each 
horizontal row corresponds to a specific rank indicated in 
the figure. This classification of the solution according to 
the ranks is comparable with the results presented in [3] . 
One can see also that the loops form "columns". Each 
column corresponds to different value of nmin which is 
equal to the number n of the number state indicated at 
the upper limit of the right border of the column. 



B. Region II: States of maximum purity 

In Fig. lUb) we have observed that the states which 
minimize the purity do not realize the whole wanted 
bound, meaning that not for all pairs of values of {pc, ^A 
we have obtained a bound on Tr {ppo). This fact is not 
surprising since the trial solution, composed by mixtures 
of number states, can achieve purity p = \, only for the 
number states. The rest of the bound which we call re- 
gion II, must be realized by states which maximize the 
purity, if pc and Tr (ppa) are given. 



In order to obtain these states we apply first the La- 
grange multipliers method for a generic density matrix 
p. Thus we obtain the extremum, i.e. minimum or max- 
imum value of the purity for a state p with given pQ and 
Tr (ppg)- In addition we require, without loss of gener- 
ality, that covariance matrix of p is symmetric in x and 
p. We express these constrains on the solution in the 
following way, 

1. The state is normalized Tr (p) = 1. 

2. The reference Gaussian state is a thermal non- 
displaced one, 

Tr (px) = Tr {pp) = Tr (pxp) = (27) 

and its purity is fixed 

Tr(p(2n + 1)) = — . (28) 
MG 

3. Fixed overlap with pc 

Tx{ppg) = Tr f p— j 

where e^ = tt-^ and N = .^.^^ the normalization 
factor. 

The solution to this extremization problem is of the 
form 

Pex = ctQ + a\x -f a2P + a^xp -|- ct/^n -|- 0156^". (29) 

However the conditions given by Eq. (P7| dictate that 
ai = a2 = as = and therefore that the solution must 
be diagonal in the eigenbasis of the harmonic oscillator. 
In other words, the extremization problem on a generic 
density matrix is reduced to the problem we have solved 
in the Sec. lIIIR and it provides the region I of the bound. 
Therefore, the method of Lagrange multipliers does not 
add anything new and in order to find the states in the 
region II maximizing purity we have to consider a pos- 
sibility of the existence of a degenerate (invariant) sub- 
space of solutions. 
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Tr(ppG) 





0.0 0.2 0.4 0.6 0.8 



FIG. 1: The region I of the bound projected (a) on the 
{/iG, Tr (ppg)} plane and (b) on the plane. The 

dotted (straight) line stands for the Gaussian states, the bold 
line for the solution of rank 2 Ea. (|26|l and the dashed line 
for the purity-bounded uncertainty relation. In part (b), the 
solid lines split the bound in vertical segments according to 
different nmin and in horzontal segments according to the rank 
fimax — ?^min of the density matrices of the solution Ea. (|20[) . In 
the figure only the ranks 2-6 and rimin : - 3 are depicted. 



1. Pure states 

Let us now proceed in a different way and try to pro- 
cess the extremization problem in a modified way to re- 
veal the degeneracy. Instead of working with the mixing 
coefficients of the density matrix, we consider them fixed 
and we make variation on the wavevectors, a mixed state 
can be decomposed on. 

To clarify the procedure we start with the special case 
of pure states where p = (-01 and — X^n"'/'" I"-) 
in an eigenbasis which for convenience we choose to be 
the one of the quantum harmonic oscillator. If we re- 
strict ourselves to states which possess as a non-displaced 
reference Gaussian state then the non-vanishing ampli- 
tudes Tpn in the superposition should be at least sepa- 
rated by three vanishing ones, i.e. \ip) = X^mV'm |*m) 
where im+i — im > 3. Since, the purity is fixed to 1 we 



have to extremize another quantity which we choose to 
be the overlap 



Tr {ppc 



e--|V) /N 



1 - MG 



(30) 



assuming in addition that the pure states are normalized 
{ijjl ^p) = 1 and of fixed covariance matrix 



l/pG^{^\{2n+l)\ij) 

= V|V'„|'(2n+l). 



The functional to variate is now the following 



(31) 



fim) = {^\ 



w 



+ ai {^\ ^)+a2 (^/'|(2n-hl) 



(32) 



After differentiating over the amplitudes ipn which with- 
out loss of generality we can consider as real numbers, 
we derive the following condition on the solution j-^) to 
the extremization problem 



(e^" + CI+ C2h) \^) = 0. 



(33) 



In this form one can see that apart from the obvious 
solution where \ip) = \n) there is one more possibility 
that was not revealed when we applied the extremization 
procedure on the density matrices, i.e., a superposition 
of two number states 



IV^) = iJm \n) + i>n+i \n-\-i) 



(34) 



with i ^ 3 (so that the covariance matrix is symmetric) 

and |-0n|^ + \i>n+i\^ = 1- 

It remains to identify the integer numbers n and i, 
Eq. ([M)) . which achieve the lowest overlap. Eg. ([50)1 . for 
a given value pc, Eq. ([3T|) . According to our search the 
states which achieve minimum overlap almost everywhere 
are of the form ipn\n) + tpn+3 \n + 3). More specifically, 
one finds that in the segment 



1 



2n-t- 1 



< 



< 



1 



2n-|-3' 



(35) 



bridging the number states |7i) and \n + 1) (see Fig[TlD), 
the states which achieve minimum are 



for 



Wa^^n\n)+^Pn+3\n + S) (36) 

1 

< Pg <rn 



2n+ 1 



and 



|?/')fc = V'«-2|r^-2)+V„+i|n+l) (37) 

for rn< PG< „ \ o 
2n -I- 3 
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where the parameter r„ is the real root of the equation 

- 2(1 + n)x^ - - 6(1 + n)a; + 3 = 

in the interval [ „''",, ■ „"'", , 1. We note here that such a 
root always exists though we do not present its form ex- 
plicitly here. 

The overlap that corresponds to the states and 
|V')f, can be easily evaluated with the help of Eqs. ([5D)) - 
(EID, 



Tr(/5/3G)a = 
and 



1 



3(1 + Mg) 
(1 - - 2n/iG) 



(7^G + 2n/iG - 1 
1 - A^G 



1 - MG 
,1 + Mg 



1 + MG 



(38) 



3(1 + /xg) 
+ (-l + 3AiG + 2n^G) 



,1 + MG 



However the state Ea. (l57|) . is not defined when 

n = or 1. In this case the minimum Tr (ppc) is rather 
obtained by the states 

IV')^ = tpn \n) + \n + 1) (40) 

with n = 0, 1 respectively. The states in Ea. (l40|) do 
not possess a un-shifted thermal reference Gaussian state 
and, apart from numerical evidence, we do not have a 
formal proof for their extremal properties. Using the 
Wigner function for the states Eq. (|40l) expressed in terms 
of Laguerre polynomials as 

Wn{x,p) = ^ U 



-X —p 



X ((|^„|'l„ (2x2 + 2p2) 

' : Re (V-nC+i - ip)) (2x2 + 2p^) 



+ 1 

we have calculated the corresponding overlaps 

Tr(/5/5G)^,o = -e'°'~^°+= 

-2a5 + Sq" - 12^3 + 5a2 + 4q; - 4 



(41) 



(2 - a)3/2 (2a2 - 3a + 2) 



5/2 



(42) 



Tr(ppG)£ 



2e''° 



2(a-l)o 



(3 

X (64ai" 
+ 6004a^ 
- 1908q2 



a)5/2 (4a- 



n9/2 



2 -5a + 3)" 
560a^ + 2156a* - 4668a^ 



42110=^ + 4940^ 
837q- 162). 



1938a-' 



(43) 



The parameter a = = 1 — IVA: 
corresponding equations for fiQ 




is defined by the 



PG,,o = -1/V2^^V4(1 - Q)a + 2a - 3 (44) 
PGf,,i = -l/\/2a - 5^8(1 - a)a + 2a - 5 (45) 



We unify the resuhs, Eqs.(I2Sl)-(|3ni), (|121)-(|1S]), and we 
present them graphically in Fig. ??. This graph sim- 
ply represents the minimum value of the uncertainty 

\J 711722 — l7i2|^ ~ 1//^G of a pure state given the 
value of its overlap Tr (ppa) (or the non-Gaussianity S). 
The line in the inset formally stands for the function 
Fip^l,6) 



711722 -I712I >F{1,6) 



(see Eq.®), which we have not explicitely derive since 
this requires an inversion of Eos. l\38\i - l\S9]\ . Interestingly 
this extended version of the Schrodinger-Robertson un- 
certainty relation is saturated not only by the ground state 
of the harmonic oscillator but by all number states. 



7pp I 7xp I 



12 


-|r«Pp 




1 


10 




1 

r 












Ci 


O.L 0.3 


0.3 11.4 





0.0 0.2 0.4 0.6 0.8 1 



:„Tr[ppe] 



FIG. 2: An extended version of the Scrodinger-Robertson re- 
lation for pure states Ea. (|III B l|l where the minimum value 

on the uncertainty '\J 711722 — |7i2|^ depends on the quantity 
Tr{pGp) (or 5). The set of minimizing states is comprised by 
the states Eg. pSI) solid blue line, \4>)i^ Ea. (|39p solid gray 
line, and Eas. (|42p - (|45p dotted line. The number states 
are included in the set of minimizing states and are marked 
on the figure as dots. 



2. Mixed states 

Moving away from the plane of pure states (density 
matrices of rank 1) it is logical to proceed with the search 
for extremal states starting from density matrices of low 
rank. The first case to be considered is that of rank 
2. According to the spectral theorem, a unique pair of 
mutually orthogonal pure states , \'ip2) always exists 
such that p = pi\ipi) {ipil + p2 IV'2) (V'2|- We consider 
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now that the probability mixing coefScients pi and p2 
are fixed and we look for the states for which the 

overlap becomes minimum. As before we assume that p 
possess a reference Gaussian state with no angular depen- 
dence on the phase space and we impose the constrains 
of normalization {ipi \ipi) — 1 and fixed covariance matrix 
Eg. (1^51) . As we did in the special case of pure states, we 
decompose the states in the eigenbasis of the har- 
monic oscillator lipi) = ^V'i.n l'^-) and apply the method 

n 

of Lagrange multipliers by differentiating over the ampli- 
tudes V'i,n- The conditions on the vectors {ipi) that we 
obtain are the following 



just present the results below, 



(e^" + C3 + C2h) \ij2) = 0. 



(46) 
(47) 



Prom the Eqs. (P5)) - ((i71) we conclude that in addition to 
simple mixtures of number states (|'0i) = K) and \i})2) = 
which we have discussed already in Sec. IIIIA. there 
is an additional solution of the form 



IV'i) = A.i N> + A,j |j) 
1^2) = |fc) 



(48) 
(49) 



where j ^ i + 3 and k ^ i, j. With simple observation we 
have arrived to the conclusion that most probably among 
the rank 2 states of this form the ones which achieve the 
minimum overlap (with some exception that we discuss 
below) are the states, Eas. p51) - P5|) . where j = i + 3 and 
k = i + loTi + 2. More precisely, the states of the 
minimum overlap are 



Pi 



■ p {{Tpn \n) + tl)„ 

{l-p)\n + i){ 



3))(Va 



31)) 
(50) 



where z = 1 or 2 and \^n\ + iV'n+al = 1- When p = 0, 
Eg. ([SO)) , the solution Eq. ([5S)) for pure states is recovered, 
while for ipn = and i = 2 (or ipn+3 — and i — 1) 
one arrives to the rank 2 states Eq. (pS|) of the Region 
I. As for pure states, for mixed states the Region II of 
the bound is not covered totally by the states suggested 
by the optimization problem but there is part that the 
(shifted) states 

P3=a \n) {n\ + {l-a)\n+ 1) {n + 1| 

+b\n) {n + l\+b*\n + l) {n\, (51) 



seem to intervene, where |6| G —y^a{l — a), y^a{l ~ a) 

and n = 0, 1. In the limiting case where \b\ = y^a{l — a) 
the corresponding states |?/')^ Eq. (|40l) are recovered from 
the Eq.([5T|). 

For the states pi, p2 the quantities of interest, Tt{pqp) 
and pg, are given by Eas. (P(I)) - (PT|) . For the states p^ the 
calculations are more involved, and for this reason we 



Tr(pG/5)3,o = c+^iM^-^ 



(2{a ~ 1) \b\^ + 2(a - 2)a \bf + (a - 2)^) 



(2-a)3/2 (-a-2|&|V2 



5/2 



(52) 



Tr(/5G/o)3,; 



ga+4|tH-3 



(3-a)5/2 (-a-4|6|V3) 
2 (^16 (4a2 - 15a + 15) \bf 
+ 2(a - 3)2 (lOa^ - 9a - 25) \bf 
-F2(a-3)^ (a^ -I- 3a- 10) \bf 
(a - 3)^(a- 2) 

-H8 (Sa^ - 45a2 + 70a - 21) 



9/2 



(53) 



where the parameters a and b are the same as in Eq. (I5ip . 
The latter connect indirectly (in a parametric way) the 
overlap with the other quantities of interest. 



/I = + (1 - a)2 + 2 



(54) 



PG3,„ = 1/^(20 - 2n - 3) (2a + 4 \bf (n + 1) - 2n - 3^ 

(55) 

with n = or 1. 




FIG. 3: The region II of the bound projected on the plane 
{/i, /ic}. We attribute to the bound three different colors; 
Blue for the part realized by the states po Eq. (|50|) . Green for 
the satets p\ Eq. (|50|) . and Red for the states pz Eg. HSip . 

If the analogous extremization procedure is applied to 
mixed states of higher rank, then the number of the ob- 
tained conditions on the solution is increasing linearly 
with the rank. For instance, for the rank 3 states which 
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can always be expressed as p ■ 
P3 IV'3) ('/'3I with (ip'] ipi) = 



= Pi IV'i) (^i|+P2 \ip2) {ip2\ + 
5ij and J^Pi ~ 1' '^'^^ 



tains in addition to 
condition 



the conditions Eqs.(|46])-(l47l), the 



For these areas clearly the function F (/i, S) in Eq.® can- 
not be defined. This problem does not appear when one 
uses the quantity Tr (ppc) instead. 



(e"''" + C4 + C2n) IV'3) = 0. 

This suggests an extreme solution of a form, e.g., pi + 
|m) (to I , where pi is defined in Eq. (j50p and \m) a Fock 
state with m ^ n, n + 3, n + i. After investigation of dif- 
ferent possibilities, we saw that such states do not achieve 
overlap lower than the rank 2 ones. This fact led us to 
the conclusion that the identified rank 2 density matri- 
ces, Eqs. ((50)) - ([5T|) are the extreme solutions which realize 
the Region II of the bound. 

In Fig. lIIIB^ in analogy with the FiglT](b), we present 
how the Region II of the bound, projected on the plane 
{/^, pg}^ is separated to three different parts, realized by 
the states pi, p2 Ea. ([50|) and ps Ea. (|5T|) respectively. 



The total bound 




FIG. 4: The non-Gaussianity uncertainty relation as a bound 
in the {pG,p,S} parametric space. The part of the bound 
that is to the right of the curly black line represents the Region 
II of the bound and the rest the Region I. In the inset the 
intersection of the bound with the plane p = 1 is presented. 

The states which minimize the purity (Region I) to- 
gether with the states which maximize the purity {Re- 
gion II) form a surface which bounds from below all the 
mixed states in the parametric space {/i, pa, Tr (ppc)} or 
equivalently bounds from above in the space {p, pc^}- 
In Fig. |4] we present the total bound that we name the 
non-Gaussianity bounded uncertainty relation. Clearly 
the minimal uncertainty, F{p,d), depends on the degree 
of non-Gaussianity of a state. Although the dependence 
is "smooth" for the low-purity region (up to purity 1/2) 
the structure of the bound becomes more complex for 
higher values of purity and for some very small regions 
the bound projected on the plane {p, 6} is double valued. 



IV. THE HUDSON'S THEOREM AND 
UNCERTAINTY RELATION FOR MIXED 
STATES 



According to the arguments in W\ the Hudson's the- 
orem for mixed states can be formulated as an upper 
bound on the non-Gaussianity for the states with a posi- 
tive Wigner function. The idea behind this formulation is 
that since among the pure states only the Gaussian states 
can possess a positive Wigner function, for mixed states 
one expects that the maximum degree of non-Gaussianity 
of a state with a positive Wigner function is lower than 
that of a generic state. In practice, this bound can be 
utilized as a set of experimentally verifiable, necessary 
conditions for a state to have a positive Wigner func- 
tion or alternatively, as a set of sufficient conditions for 
a state to posses a non-positive Wigner function. Such 
conditions can be essential in the construction of "non- 
classical" states which are^useful resources in protocols 
of quantum information [14]. 

Even though it is not trivial to identify the whole 
bound in the space {p, pc, 5} for states with non- negative 
Wigner functions, we can conclude using the same argu- 
ments as in Sec. IIIIt \ that the purity of the states with 
non- negative Wigner function for given pQ and Tr (ppg) 
is minimized by mixtures of number states. To our 
knowledge there is no necessary and sufficient condition 
that one can impose on the coefficients of a mixture of 
number states to guarantee positivity of the Wigner func- 
tion. Due to this difficulty we cannot proceed with fur- 
ther analytical identification of the states which realize 
the Hudson's theorem for mixed states. 

On the other hand from the results in Sec. IIIII we can 
numerically obtain some information about the Hudson's 
theorem for mixed states. To do so, in the bound FigH] 
for all states we identify numerically the part realized 
by states with positive Wigner function. We depict this 
common area, where the bounds for all states and for 
states with strictly positive Wigner function are overlap- 
ping, as a shaded region in Fig|S] In ^TE] we have gone one 
step further and we have identified numerically the whole 
region I of the bound for the states with positive Wigner 
function. We have concluded that the two bounds stay 
very close when the purity of the state is low (for a fixed 
uncertainty) and we obtained an indication that the dis- 
tance between the two bounds is increasing with purity. 
In addition in [15| the lowest degree of Tr (ppc) that a 
state with strictly positive Wigner function may achieve 
for a given uncertainty, has been identified. 
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6 




FIG. 5: The part of the bound (red) reahzed by states with 
positive Wigner function projected on the planes (a) {no, 5} 
and (b) {/xc/i}. In figure (a) it is plotted together with 
the upper non-Gaussianity bounds for all mixed states (solid 
line) achieved by pure states and the boundary line (dashed 
line) among the Region I and // realized by the rank 2 states 
Eq.((26|. 



V. CONCLUSIONS 

Inspired by the discussion in [3] and some recent re- 
sults on the extension of Hudson's theorem for mixed 
states 0, we have extended the purity-bounded uncer- 
tainty relation, adding one more parameter, namely, the 
degree of non-Gaussianity. According to the results pre- 
sented here, the minimal uncertainty for quantum states 
strongly depends not only on the degree of mixedeness 
but also on its non-Gaussian character. The set of 
states which saturate this uncertainty relation includes 
the Gaussian states. From this we can conclude that 
even though Gaussianity is not a requirement for the 
minimization of this uncertainty relation, Gaussian states 
are still included in the set of minimizing states and 
they dominate as purity tends to one. Interestingly, the 
non-Gaussianity bounded uncertainty relation is valid for 
pure states and there it is saturated, among others, by 
all eigenstates of the harmonic oscillator which appear as 
extremal points of the bound. 

Moreover, by identifying the states which possess a 
positive Wigner function among the ones which satu- 
rate the non-Gaussianity bounded uncertainty relation, 



we come closer to the extension of Hudson's theorem for 
mixed states. Gaussian mixed states (thermal states) 
both saturate the non-Gaussianity uncertainty relation 
and satisfy the extended Hudson's theorem. We should 
note that in we employ numerical methods to compli- 
ment the results presented here concerning the Hudson's 
theorem. The main result in [Tsll is the identification of 
the minimum degree of trace overlap that a mixed state 
of fixed uncertainty may achieve. 

In this work we have chosen purity as a measure of the 
degree of mixedeness of the state. If the von Neumann 
entropy is employed instead [8] , then the results are sim- 
ilar to the one presented in this work and one arrives to 
the non-Gaussianity extension of the "entropy-bounded" 
uncertainty relation suggested by Bastiaans f&l . The ad- 
vantage in using the entropy, is that the solution is au- 
tomatically a positive matrix and one does not need to 
employ the ansatz we have used in the current work, to 
impose positivity (Ea. (jl9p ). Furthermore, the solution 
to the extremization problem contains two "branches" , 
one bounding the trace overlap from above and one from 
below for fixed uncertainty and purity. Note than in the 
current work we derive only the lower bound on the quan- 
tity of trace overlap. On the other hand, the states which 
minimize the extended uncertainty relation in Q are 
mixtures of all number states with n — > oo, the solution 
cannot be expressed in a closed form and one needs the 
help of numerical methods to visualize the results. 
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Appendix A: A parametric expression for the 
purity-bounded uncertainty relation 

Here, we derive an expression for the purity-bounded 
uncertainty relation alternative to the one presented in 
Ref. [5]. The derivation of the non-Gaussianity bounded 
uncertainty relation and its approximate formula in 
Sec. mil is a direct generalization to higher dimensions 
of the procedure we present here. 

In [5] it is proven that states with angular-independent 
Wigner function minimize the relation. We use this result 
and we start with states which are finite convex combi- 
nations of the number states 

N 

p=Y,Pn\n}{n\ (56) 

n=0 
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where 



N 

E 

n=0 



Pn = 1 



(57) 



and < P„ < 1 • The purity of p, and the purity of the 
reference Gaussian state are determined by 



N 



n=0 



N 

E 

n=0 



Pn{n + 1/2) = l/2^c 



(58) 
(59) 



respectively. As in j4| we use the Lagrange muhiphers 
method and we arrive to the foUowing form for the co- 
efhcients of the states which minimize the purity fi[p\, 



(60) 



where Aiand A2 are to be determined by the conditions 
Eas. ([571) and (|59p . However the parameter N, the upper 
hmit in the summation, remains undetermined and to re- 
solve this ambiguity we introduce a continuous parameter 
y that satisfies 



Ai+A2y^0. 



(61) 



This additional equation permits us to define N as the 
integer part of y, iV = [yj . 



We employ the conditions, Eqs. ((57|) and ((6T|) , and 
express Aiand A2 in terms of N and y. Then we sub- 
stitute the result into the Eqs. ([58| . ([59l) and we obtain a 
parametric relation for the extremal solution provided by 
the Lagrange multiplier's method, 



3 {N - 2y) 



N{5- 
2 (N- 



47V) - 
-2N^ 



6(U 

- 6Ny 



N)y' 



3{l + N) {N-2yY 



(62) 
(63) 



where N = [y\ and y G [1, 00). We note that by setting 
N = y one can arrive to an approximate formula 



9 - (Mg=) 



that reproduces the curve of Eas. ([5^ - (I55|) in very good 
approximation when < /ig < 3/5. For the rest of 
the region, 3/5 < fia < 1, one should better employ the 
exact solution 



2y.l 



(64) 



The obtained results Eas. ([5^ - (p5|) give the lower 
bound (see Fig. 1(a)) on the plane jic and /x and are 
equivalent to the expressions derived in 



Appendix B: Parametric expression for the non-Gaussianity Bounded uncertainty relation 



Here we present the explicit expressions for the trace overlap and the purity that realize the part of the non- 
Gaussianity bounded uncertainty relation covered by mixtures of number states. 



pT =- {^—^^ f 2Air + 2^2n,„i„r - A2r + A2 + A^r'' f 

\VTW (^) ^ ^ (' (^'^' " (^) ^ ^ - 2Ai(r + 1)^ - A2{r + l){2lx\ r + r + 1)^ 

At"'' = - ^ (nniin - [x\ - 1) (6^1^ + QA1A2 (n„,i„ + [x\) + (2n^in [x\ + n„iin(2n„i„ - 1) -t- 2 [xj V ) ) 

(65) 

- ^37^;—^ (^] ^ ^ (Mr - l)r (^] ^ ^ - 4Ai(r + 1)^ - 2A2{r + l)(2LxJ r + + 1) 



Ar{r + l)^ \r + I J \ \ \r + I 



(66) 
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where r stands for /ic x for X2 and Ai, A2, and are defined by the hnear system of equations 

Ai + A2X + A^r (^^^^ /(r + l)=0 (67) 

Ai(-n„,i„+ [a;J +1)+A2((1 -n,ni„K„i„+ N(W +l))/2 + A3 ||(^^^^ " (^Tt) ) = ^ (68) 
^1 (-n,„i„2 + N V 2[a;J + l) + ^^2 (-4n„,in^ + Sn^in^ + n,ni„ + lx\ (4[xJ V 9La;J + 5 

+ A3((^] ""'°(2n„i„r + l)- fi-^V (2(La;J + l)r + 1) ) /2r = 1/r. (69) 



I 1 / V mill- ' I 1 1-1 

r + ly \7'+l 

One should first solve the above system of linear equations to express A's in terms of the parameters r, x and n^in and 
substitute the resulting expressions in Eqs. (l65l) - (|66p . One should then attribute an integer non-negative value to rimin- 
In this way the running parameters left are r G [0, 1] and x e [2, 00]. However their limits are furhter constrained by 

/ Y J, \ "mill 

Ax + vlarimin + A^r [ \ / (r + 1) > 0, V (70) 

1 — r ^ 



Ai+A2(n^in-1) + A3r ( ) /(r + l)<0, z/ ^ (71) 



which ensure that n,„i,, remains constant 
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